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Abstract 

The lattice dynamics of solid 4 He has been explored using pulsed NMR methods to study the motion of 
3 He impurities in the temperature range where experiments have revealed anomalies attributed to superflow 
or unexpected viscoelastic properties of the solid 4 He lattice. We report the results of measurements of the 
nuclear spin-lattice and spin-spin relaxation times that measure the fluctuation spectrum at high and low 
frequencies, respectively, of the 3 He motion that results from quantum tunneling in the 4 He matrix. The 
measurements were made for 3 He concentrations 16 < x% < 2000 ppm. For 3 He concentrations X3 = 16 
ppm and 24 ppm, large changes are observed for both the spin-lattice relaxation time T\ and the spin-spin 
relaxation time T2 at temperatures close to those for which the anomalies are observed in measurements of 
torsional oscillator responses and the shear modulus. These changes in the NMR relaxation rates were not 
observed for higher 3 He concentrations. 
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L INTRODUCTION 

The observation of distinctive anomalies in the frequency and dissipation of torsional oscillators 
containing solid 4 He, often referred to as non-classical rotational inertia fractions (NCRIFs), by 
Kim and Chan 2 3 has stimulated enormous activity. This is because the NCRIFs could be the sig- 
nature of a supersolid state as outlined by Leggett 4 and predicted by Andreev and Lifshitz several 
years ago. 5 A large number of independent experiments 6-10 have shown that the NCRIF magnitude 
and temperature dependence are strongly dependent on defects such as 3 He impurities 1112 and the 
quality of the crystals, and can be made very small by very careful annealing. 3 The observation of 
NCRIFs was preceded by reports of a striking anomaly in the measurement of sound attenuation 
in solid 4 He by Goodkind and colleagues who had compared their observations to that predicted 
for a true phase transition of a fraction of the solid to a supersolid state. 1314 Additional support 
for the interpretation of the NCRIF anomalies in terms of a true thermodynamic phase change to a 
new state of matter was provided by the observation of a small but distinct contribution to the heat 
capacity of the solid at the same temperatures as those reported for the onset of the NCRIFs. 15 

A straightforward interpretation in terms of a transition of part of the solid to a coherent su- 
perfluid component has, however, been hampered by a number of puzzling observations. These 
include the lack of evidence for a critical exponent, an apparent very low critical velocity 7 , the 
absence of experimental evidence for fourth sound modes 1617 and the null results of attempts to 
observe pressure- induced superflow though small restrictions, 18 although mass transport has been 
observed with the use of porous vycor glass conduits used to allow mass flow to the solid via su- 
perfluid in the pores. 19 ' 20 In addition, measurements of the shear modulus of solid 4 He by Beamish 
and colleagues 2122 have revealed a prominent frequency dependent change in the elastic shear 
modulus with an enhanced dissipation peak having a temperature dependence comparable to that 
observed for the NCRIF. These results suggest that the torsional oscillator anomalies may result 
from unusual elastic 23 or viscoelastic 24 properties of the solid 4 He rather than superfluidity. Other 
interpretations attribute the anomalies to macroscopic superflow mediated by defects or disloca- 
tion networks 25 or in terms of a vortex model 26 for which the high temperature tail of the NCRIF 
is associated with the finite response time of vortices to the oscillating flow fields in the TOs. An- 
other set of researchers suggest that the superfluidity arises from non-equilibrium behavior leading 
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to superflow along defects or the formation of a quantum "superglass" around extended defects 
with ultra-slow relaxation dynamics reminiscent of glass dynamics. 28 

The need for any interpretation to account for both sample dependent coherent non-inertial flow 
and changes in the lattice elastic properties was confirmed strikingly by Kim et al. 29 who observed 
that the TO response could be changed significantly by rotating the cryostat, with the resonant 
frequency changing with speed, but that at the same time, the TO mode that showed significant 
drive dependence was not susceptible to changes in the elastic modulus of the lattice. 29 

It is clear that the dynamics of the 4 He lattice plays an important role in the low temperature 
bulk properties of solid 4 He and rather than observing a phase transition to a supersolid state one 
may be observing a complex thermally excited dynamical response, or both, in overlapping tem- 
perature ranges. It is therefore important to study the microscopic dynamics of solid 4 He using 
different techniques and especially non-invasive techniques since all previous studies have in- 
volved the application of macroscopic external mechanical stimuli. Measurements of the quantum 
tunneling of the 3 He atoms in the solid provide a unique means of probing microscopic lattice 
properties because the 3 He - 4 He exchange rates depend exponentially on the lattice separation and 
on the magnitude and dynamics of the crystal field deformation that surround impurities. 30 ' 31 The 
characteristic NMR relaxation rates are determined by the modulation of the nuclear dipole-dipole 
interactions by the tunneling motion and the scattering of the diffusing atoms by the crystal de- 
formation field around the 3 He impurities and other lattice defects. 3132 The NMR relaxation rates 
are therefore very sensitive to the local elastic properties of solid 4 He and to any changes in the 
crystal ground state that would modify the tunneling rate. We have therefore carried out systematic 
measurements of the nuclear spin relaxation rates of 3 He impurities in solid 4 He for a wide range 
of concentrations with an emphasis on low concentrations x 3 < 30 ppm for which one expects 
well characterized NCRIFs from previous studies and yet sufficient to obtain good signal to noise 
ratios for NMR measurements. Also for these concentrations, the amount of the 3 He localized 
on dislocations estimated from the expected density of dislocations for well-annealed samples 33 
is expected to be a small fraction of the total, less than 1 ppm, so that one can be sure that the 
experiments probe the motion of the 3 He atoms through bulk solid 4 He. Preliminary results of 
these NMR studies have been reported elsewhere. 34 ' 35 
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II. EXPERIMENTAL ARRANGEMENT 



The NMR cell was designed as a nested cross-coil arrangement with receiving and transmit- 
ting coils orthogonal to each other and the applied magnetic field (see Fig. 1 ). The inner coil is 
a cylindrical receiving coil wound around a polycarbonate sample cell which contained an in situ 
Straty-Adams pressure gauge. The outer coil is a thermally isolated orthogonal transmitting coil 
for the RF pulses that slides onto the receiving coil. This nested arrangement 36 provides (i) mini- 
mization of the unwanted pick-up of the transmitter pulse by the receiving coil, and (ii) adequate 
thermal isolation from the heat generated in the transmitting coil that is heat sunk to the still of the 
dilution refrigerator. Thermalization of the sample was assured by thermal contact with a silver 
post extending from the dilution refrigerator and the temperature was measured using a Cernox 
resistance thermometer calibrated against a 3 He melting curve thermometer. The sample gas was 
admitted via a capillary sealed with epoxy (stycast 2585 GT) at the opposite end of the cell. The 
sample cell could withstand pressures up to 100 bar and remain superfluid leak-tight. 

Samples were grown by the blocked capillary method. Gas samples of pre-determined 3 He 
fraction were prepared at room temperature by mixing pure (99.99%) 3 He and 4 He, compressing 
the gas mixture to about 50 bar and filling the sample cell with the gas mixture to a pressure of 
46 bar, after which the cell was cooled to 1 .2 K. When the temperature is cooled to the melting 
curve, the helium pressure follows the melting curve while the mixture solidifies, after which the 
pressure measured by the gauge in the cell remains constant indicating that the sample is solid. 
The samples were prepared to all have the same final pressure of 27.75 ± 0.05 bar at 1.2 K. The 
molar volume of the samples were determined to be V m = 20.85 ± 0.05 cm 3 from the PVT data 
of Grilly and Mills 37 and the formula derived by Mullin. 38 Finally the samples were annealed for 
24 hours just below the melting point except for one case for which a minimum annealing of 0.5 
hours was carried out to determine the effect of the crystal quality on the nuclear spin relaxation 
times. 

The measurements of the nuclear spin-spin and spin-lattice relaxation times were carried out 
using a superheterodyne pulsed NMR spectrometer operating at a Larmor frequency of 2.05MHz. 
At this frequency the calculated relaxation times due to exchange motion were of the order of 
10 4 s. at 3 He concentrations x 3 = 20 ppm, extrapolating from previous measurements. 39 ^ 2 The 
longitudinal relaxation time is expected to increase rapidly with frequency 39 , making studies at 
higher Larmor frequencies extremely difficult. 
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The signal/noise ratio for samples with X3 ~ 20 ppm is very weak (< 1CT 2 ) for a 10 kHz 
bandwidth employing a standard geometry with resonant NMR circuit connected directly to a 
room temperature amplifier. In order to enhance the signal/noise we developed a preamplifier 
that could be operated at low temperatures in the applied magnetic field. 36 The device used a 
pseudomorphic high electron mobility field effect transistor with adjustable bias so that the output 
could be matched to a 50 ohm cable. The total power dissipation was 0.5 mW and to attain the 
lowest sample temperatures the amplifier was separated from the sample cell by a 75 cm length 
of cable and anchored to a IK thermal plate. A typical NMR echo recorded using a 90 x — I80 x 
RF pulse sequence for a 500 ppm 3 He sample is shown in Fig. 2 after signal averaging for 10 
pulse sequences using a 17 kHz bandwidth. As it is necessary to wait several times Ti between 
pulse sequences, up to 15 hours were required to obtain each data point at the lowest temperatures. 
Therefore the total number of data points is modest, but sufficient we believe to reveal and confirm 
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FIG. 1. (Color online) Schematic representation of the low temperature NMR cell. The sample in this 
NMR cell is cylindrical, 5.0 mm diameter and 7.5 mm long. The preamplifier and tuning capacitor are 
located on a 1 K cold plate located a distance of 1.2 m from the sample cell. The RF transmitting and 
receiving coils are simplified in this figure. 
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relevant features in the temperature dependence of the relaxation times. 




Time (ms) 



FIG. 2. (Color online) Example of NMR Hahn echo for a sample of 500 ppm 3 He in solid 4 He at T=0.4K, 
showing the in-phase and out-of-phase components. 

One of the most important properties of the sample that must be determined is the exact 3 He 
concentration, because for the long capillary line used (1.2 m below the IK plate) an appreciable 
fraction of the 3 He atoms can plate out on the capillary walls and lead to changes in the 3 He 
concentration of the solid sample compared to that of the original gas concentration. Fortunately, 
the NMR signal itself allows one to obtain an accurate calibration by measuring the NMR echo 
amplitude for fixed spectrometer gain at intermediate temperatures and comparing to a relatively 
high concentration sample (2000 ppm) for which any fractional change would be small. In Fig. 3 
we show the temperature dependence of NMR echo amplitudes for several different samples. 
All the samples show the characteristic Curie law behavior until reaching temperatures of the 
order of 0.1K where phase separation occurs. At these low temperatures the 3 He forms Fermi 
liquid droplets with temperature independent NMR amplitudes and with different relaxation times. 
Using this procedure one can measure the 3 He concentration to approximately 5% accuracy. For 
example, for a sample prepared from a nominally 30 ppm gas mixture, the final solid 3 He fraction 
was measured to be 24 ppm. Note that a small deviation was observed in the amplitude of the echo 
signal for the 16 ppm sample at T = 175 mK and this is attributed to the sudden and unexpected 
sharp increase in T\ at that temperature. Waiting times between measurements were set at 5 times 
T\ to provide 5% accuracy in amplitude measurements but at 175mk this drops to 2.5 times T t 
and a small decrease in observed amplitude results. 
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FIG. 3. (Color online) Observed temperature dependence of NMR echo amplitudes for several different 
samples measured for the same spectrometer pulse and gain settings. The broken lines represents the Curie- 
law behavior. 

III. EXPERIMENTAL RESULTS 




10 100 1000 

Temperarture (mK) 



FIG. 4. (Color online) Observed temperature dependence of the spin-lattice and spin-spin relaxation time 
for two samples (i) X3=16 ppm, and (ii) X3=500 ppm. 



Fig. 4 shows the observed temperature dependence of the nuclear spin-lattice relaxation time 
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(Ti) and the spin-spin relaxation time (T 2 ) for a sample with x 3 = 16 ppm and the comparison 
with the results for a sample with x 3 = 500 ppm. The temperature independent tunneling induced 
relaxation occurs for 0.2 < T < 0.55 K and a small upturn in Ti is observed for T > 0.55 K with 
a corresponding downturn in T 2 . This high temperature behavior is attributed to the onset of the 
effect of thermally activated vacancies. The most interesting features are the sharp peak in Ti at 
T = 175 mK and a less well-defined minimum in T 2 at the same temperature for the sample with 
£3=16 ppm. The same behavior as a function of temperature is observed for a sample with x 3 =24 
ppm prepared with a considerably reduced annealing time. 35 . The sharp changes in Ti and T 2 are 
not observed for samples with concentrations x 3 > 200 ppm. These features near 175 mK are un- 
expected in the traditional interpretations of the relaxation of dilute 3 He in solid 4 He that attribute 
the relaxation to a temperature independent quantum tunneling J 34 of the 3 He atoms through the 
4 He lattice. The temperatures at which these sharp features for Ti and T 2 occur coincide with the 
temperatures for which anomalies are observed for measurements of the rotational inertia 2 ' 3 and 
shear modulus experiments 22 ' 43 . On cooling below 175 mK, a precipitous drop in the value of Ti 
is observed at T = 95 mK which is the temperature associated with the phase separation of the 
solid mixture into 3 He rich droplets in otherwise pure solid 4 He. This interpretation is confirmed 
by the observed change in the NMR amplitude below the phase separation temperature as shown 
in Fig. 4. and by the significant hysteresis on cycling through the phase separation temperature. 35 



IV. DISCUSSION 

The theoretical treatments of the relaxation of dilute 3 He atoms in solid 4 He 31 ' 32 44 ^ 6 all take 
into consideration the lattice deformation surrounding a 3 He impurity due to the increased zero- 
point motion of the 3 He atoms with respect to the lattice compared to the 4 He atoms 47 The resulting 
deformation has been described in terms of a long range anisotropic interaction between 3 He atoms 
given by 

K(r tj ) = K (3cos% 1 -l)(^) 3 (1) 

for the deformation energy at site j due to an impurity at site i, where ao is the lattice constant. 30 ' 48 . 
The 3 He atoms travel through the lattice and scatter from one another with a closest approach 
given by the distance b c for which the kinetic energy of the tunneling particle (~ J 34 ) becomes 
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comparable to the elastic deformation energy 



(2) 



The mean separation of the 3 He atoms r rn = a x 3 and for x 3 > £30 = ^ ~ 1CT 3 the 3 He 
atoms are in continuous interaction with elastic fields of one another. For x 3 < x 30 , the atoms 
move coherently by tunneling until scattered by other 3 He atoms or other lattice defects. This 
simple estimation of the concentration £30 for which there is a cross-over from simple coherent 
diffusion to the continuous interaction regime was analyzed in detail by Landesman and Winter 45 
using a moment expansion for the calculation of the physical diffusion and they found a very 
different value with x 30 ~ (j^) 2 ~ 1CT 6 . Huang et al. 32 re-examined this estimate using a more 
precise accounting of of the energetics of the scattering process and estimated £30 ~ 1CT 4 . As 
we will show below, the Landesman model gives a good description of the values of the spin 
lattice relaxation time measured by different research groups with a characteristic concentration 
dependence of Xi oc x 3 2 ^ 3 for x 3 > 10~ 4 and a very different concentration dependence below 10 
ppm. 

In order to understand the possible origins of the experimentally observed peak in 7\ at low 
temperatures we need to examine how the tunneling motion of the 3 He atoms and their accompa- 
nying lattice deformations determine the nuclear spin-lattice relaxation. Landesman 31 treated the 
motion of the 3 He atoms in the presence of the deformations using a fictious spin model. In this 
model a fictious spin Si takes the values or ±1 according to whether a site i is occupied by a 4 He 
atom or a 3 He atom with real nuclear spin I? = ±1, respectively. The probability that a site j is 
occupied by a 3 He atom is represented by tj = (S 3 Z ) 2 and the elastic deformation energy by the 
Hamiltonian 



where P jk = ^(S^S^S^S^, + S^SyS^S^,) + h.c is the permutation operator for atoms at sites i 
andj. 

The time dependence of the fictious spin operators Sj and Tj are given by 




(3) 



The tunneling Hamiltonian is given by 




(4) 



S+(t) = e iIiKt Sfe 




(5) 
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where Wj = 2 J2j<k Kjk T k- hi this time dependence Landesman did not include a lattice relaxation 
term to account for relaxation of the lattice as the deformation surrounding the impurity atom 
moves through the lattice. Given the recent results from elastic studies of solid helium at low 
temperatures this could be an important effect. The relaxation times can now be estimated by 
calculating the correlation functions of the nuclear dipole-dipole interactions. We consider the 
reduced correlation function 

rg« = ^(WiJfW (6) 

where the TV, are the irreducible nuclear spin operators (that transform analogously to the spherical 
harmonics K, m (%)). Tff(t) = e iHMt T m e -iH u t with H m = H K + H T the total Hamiltonian 
responsible for the motions in the lattice. For the dilute lattice the tunneling Hamiltonian does not 
commute with all and one finds that the derivative 



r m (t)/T(0) = -12x 3 J 2 u F vp (t)cos[(K tl - K jp )t}(^-) 6 (7) 



with F ip (t) = Y[ s cos(u 3 (t) where u s = x 3 (K is — K ps )t. Evaluating F{t) using the statistical 
method of Anderson 49 yields 



r(t)/r(0) = -48A x 3 Ji 4 exp[-x 3 A(K ty' 4 } (8) 

where A' = 0.28 and A = 8.77. 31 ' 45 Integrating Eq'n (8) Landesman 31 finds a correlation time t c 
given by 



-1 _ D J 34 ^-1/3 

' c 

with B = 23. 



Tc 1 = B-^x,^ (9) 
A 



The correlation time r c < u L 1 for all the experiments reported and we therefore find for the 
calculated nuclear spin-lattice relaxation time in Landesman's model 

T "=A^ 13 <l0) 

where M 2 is the second moment for pure 3 He and J e // = J^ 4 /K . As shown in Fig. 5 this 
calculated value provides a good description of the observed values of Ti in the temperature- 
independent "plateau" region for 3 He concentrations x 3 > 100 ppm for ^J- = 1.2 kHz and no 
other adjustable parameter. At lower concentrations a much stronger concentration dependence is 
observed with Ti oc x 3 . This behavior at low concentrations is attributed to the crossover from 
the continuous interaction regime to a region of coherent diffusion for which the characteristic 
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FIG. 5. (Color online) Observed 3 He concentration dependence of the nuclear spin-lattice relaxation time 
for dilute 3 He in solid 4 He. Orange circles, Schratter et al. , up triangles, Allen et al. , diamonds Schuster 
et al. 51 , down triangles 52 , open squares 53 , solid squares 35 . 

time is determined by the time for a 3 He atom to travel the mean distance r m = x^^ao between 
3 He atoms. This time is given by r coh = r m /v g where the group velocity v g = a zJ 34 . We find 
T coh = {zJ^^x^ leading to a calculated longitudinal relaxation time T^ oh = 2.8 x 10~ 3 x^ 3 
shown by the solid red line of Fig. 5 for J 34 / (2tt) = 1.2 MHz. 

The predicted temperature independent relaxation for x 3 = 16 ppm is shown by the broken lines 
in Fig. 6. While the fit to experimental data is good at high temperatures the anomalous relaxation 
observed at 175 mK is more than a factor of two larger than the Landesman prediction for a 
simple tunneling motion of the 3 He atoms through the lattice. The anomaly occurs well above the 
phase separation temperature (90 mK) and is not attributable to the phase separation which would 
lead to a sharp decrease in the relaxation time and not the sharp peak that was observed. This 
conclusion was confirmed by studying a sample with 24 ppm that was never never cooled to the 
phase separation temperature (see Fig. 6). 

It is important to note that a simple change in the elastic field interaction K associated with 
thermal population of excited states with different values for K would not explain the observed 
results and would simply lead to a broad step in the value of Ti. The results therefore imply that 
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the anomalous peak observed for T\ results either from a change in the dynamical properties of 
the lattice that occur at 175 mK or that there is a phase transition that occurs at that temperature 
and the peak observed in T\ results from the fluctuations at the phase transition. 

10000 
7500 
H 5000 
2500 


100 200 300 400 

Temperature (mK) 

FIG. 6. (Color online) Comparison of the values of the observed nuclear spin-lattice relaxation times for 
£3=16 ppm (red circles) and x%=24 ppm (blue triangles) with the theoretical model of Landesman (broken 
lines a and b) for tunneling 3 He impurities in the 4 He lattice with a fixed lattice distortion. The sharp drop 
at 90 mK marks the well known phase separation with the formation of pure 3 He nanodroplets. The 24 ppm 
sample was never cooled below 120 mK to avoid hysteresis and memory effects associated with the phase 
separation. 




A. Effect of a transition to supersolid state on the NMR relaxation times 

We now consider what would be expected if there was a simple phase transition to a state in 
which there was a small fraction of superfluid condensate, and then we will discuss how the effect 
of a dynamical relaxation of the lattice as described by Beamish and colleagues would influence 
the nuclear spin lattice relaxation time T\. 

If there was a sharp phase transition at a fixed temperature T c below which a superfluid state 
appears one would expect a sudden change in the spectral density at Tc that determines T\ and 
T 2 because of the sharp increase in fluctuations just below T c that would modulate the spin- 
spin interactions (assuming that the 3 He atoms are carried by the superfluid component). One of 
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FIG. 7. (Color online) Comparison of the observed temperature dependence of the nuclear spin-lattice 
relaxation time T\ with the dependence expected for a true phase transition at T=175 mK for a sample with 
£3 = 16 ppm as described in the text. 

the most significant results of Chan et al? is that the critical velocity attributed to the superfiuid 
interpretation of the NCRIFs is very low, typically 10~ 4 m/s. The fluctuations would therefore add 
weight to the spectral density at low frequencies and not at high frequencies (>MHz). Because 
the spectral densities are determined by the fluctuations of the 3 He dipole-dipole interactions (and 
no other magnetic interactions are present) they are normalized and thus the increase in spectral 
weight at low frequencies would result in a decrease in the spectral weight at high frequencies 
with an increase in T\ and a concommittant decrease in T 2 . 

If we designate the component of the spectral density due to 3 He particle exchange as Je(cu) 
and the component due to a possible superfiuid component as Js(w), we have 



where M 2 is the NMR second moment. For liquid 4 He the coherence length scales as £ = a (T c — 
T)~ 2 / 3 where Tc is the critical temperature, and because of the very low critical velocities (inferred 
from torsional oscillator measurements), we anticipate a characteristic relaxation time that behaves 
as t s = t so (T c — T)~ 2 / 3 for T < T c , using t S o — r m /v c , where r m = x _1 / 3 a is the mean 
separation between atoms, a is the lattice spacing and v c is the critical velocity. If the supersolid 
component component has weight a (T) then the exchange component must have weight 1 —a(T) 




(ID 
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so that the sum rule J ^- = ixM 2 is obeyed. Thus the observed relaxation rate will be given by 

7^— = a(T)J s (cj L ) + (1 - a{T))J E {u L ) (12) 

1 lObs 

where Js(wl) is negligible. If T 1E is the calculated value of 7\ for exchange motion only, we find 

T\obs 1 



T 1E ' [1 - a(T)] 



(13) 



where for a supersolid fraction we take a(T) = a x for T < T c and a(T) = for 

T > Tc. ai is an adjustable parameter related to the '"supersolid"' density which is of the order 
of 0.01 but sample dependent. 

An example of a fit using the above expression for T 1Gfes is shown by the broken line in Fig. 7 
for Tc= 175 mK and a\ = 0.01 with no other adjustable parameter. In reality we would expect to 
have a small Gaussian spread (of about 7 mK) in the critical temperature because of the inhomo- 
geneities and it would not be difficult to obtain a good fit to 7\ for reasonable numbers but no other 
experiment points to a sharp transition (except for the onset of NCRIFs in very pure samples), e.g. 
no sharp peak in the heat capacity is reported. A more stringent test of this interpretation is pro- 
vided by examining the results for the T 2 measurements for the same sample. The same argument 
as used for evaluating T\ needs to be followed for T 2 but this time we need to know the values for 
t s0 and the exchange time te- 

First we consider the theoretical predictions for T 2 for small x 3 . The values of the correlation 
time calculated by Landesman (Eqn. 10) that give a good fit to the T\ data do not lead to values of 
T 2 that are in agreement with the experimental data. This discrepancy has been reviewed by Kim 
et al. 54 who showed that the detailed spectral density at low frequencies cannot be described by a 
single Lorentzian spectral density. A better approach to evaluating T 2 was given by Huang et al. 32 
who calculated the scattering of tunneling 3 He impurities in the elastic crystal fields using detailed 
energy conservation and found T 2 H = 1.69 x 10~ 4 a;3 s. for x 3 > 100 ppm. This result is compared 
with the experimental data and Landesman's theory in Fig. 8. For very dilute concentrations, x 3 < 
100 ppm, the model of coherent diffusion as discussed above gives T 2 oh = 1.68 s. 

For these dilute samples the tunneling motion of the 3 He impurities would lead to a temper- 
ature independent relaxation, but once again, we observed strong deviations from temperature 
independence below 175 mK. However, unlike the T\ results which show a clean peak, there is 
considerable scatter in the T 2 data and the dependence observed is that of a small peak followed 
by a strong dip, and then at even lower temperature there is an order of magnitude drop at the 
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FIG. 8. (Color online) Observed 3 He concentration dependence of the nuclear spin-spin relaxation time 
for dilute 3 He in solid 4 He for several concentrations reported by different authors. The solid lines T 2 L and 
T 2 H are the predictions of the theories of Landesman 31 and Huang et al. 32 , respectively. T 2 o}l is the value 
of T2 for coherent diffusion where the mean separation is greater than the scattering length. See legend of 
Fig. 6 for definiton of symbols. 

phase separation temperature where Fermi liquid droplets form. We now calculate what would 
be expected for a simple superfiuid phase transition. In terms of the spectral densities (at zero 
frequency) for the superfiuid and normal components we have 

1 



a(T)J s (p) + [l-a(T)]J E (p) 



while the fixed temperature independent tunneling term is 

1 



(14) 



T, 



2E 



We therefore find 



T, 



20bs 



Mo) 



1 



(15) 



(16) 



T 2E ba(T) + [1 - a{T)\ 
where b = Js(0)/ Je(0) = tsq/te > 1 and t e is the characteristic tunneling time in the normal 
solid. Because of the large value of b, the result for T 2 obs does not look like the inverted image of Ti 
(as one would naively expect) but is much wider. An example of an attempted fit for the adjustable 
parameters 6 = 4 and a\ = 0.01 is shown in comparison with experimental values of T 2 in Fig. 9. 
Even allowing for the scatter of the data, this approach does not provide a good description of 



15 



the temperature dependence of T 2 . We should note that if the system is very inhomogeneous with 
different parts of the sample having different values of b and a\ one would expect a lumpy spectral 
density at low frequencies and a corresponding scatter in the results for T 2 , but since T\ measures 
the spectral weight at MHz frequencies, one could still observe a sharp peak for T\ because it 
is the integrated weights of the low frequency components of the spectral densities that must be 
subtracted from the high frequency spectral density. In this sense the T\ measurements are a much 
cleaner test for the existence of a phase transition that the T 2 measurements. 
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FIG. 9. (Color online) Comparison of the observed temperature dependence of the nuclear spin-spin 
relaxation time T 2 for a sample with X3 = 16 ppm with the dependence expected for a true phase transition at 
T=175 mK (broken purple line) as described in the text with the high temperature lmit given by Landesman's 
model. 



B. Effect of lattice dynamics on the NMR relaxation rates. 

An alternative approach to interpreting the NMR relaxation rates is to assume that as a result 
of the tunneling of the 3 He atoms there is an internal stress on the 4 He lattice as the distortion 
around the 3 He impurity moves with or attempts to move with the impurity as it tunnels from site 
to site. The mean frequency of this motion is co x = x^ 3 K w 3.5 x 10 3 Hz. We can surmise that 
the relaxation of the lattice due to this time dependent disturbance is long and adds a bottleneck 
for energy exchange between the tunneling excitations and the thermal bath of phonons. The 
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correlation function for the lattice operators r(T) of Eq'n (3) can be written as 

(r(t)r(O)} = (r(0) 2 } R exp(-iu x + r^t (17) 

where ( )r represents an average over lattice coordinates, r s is the lattice averaged relaxation 

time which we assume corresponds to the Debye relaxation r s = r exp(E /T) used to interpret 
the results of the shear modulus measurements. 2123 Integrating Eq'n (17, we find for the additional 
lattice relaxation time 

Ta; ( T ) = ri TT~( eJt^ = — — 2 (18) 

1 + (UxToe^ ' 1 y u x l + u 2 

where u = u x r e^ E °^ T \ r and E are adjustable parameters that we expect to be comparable to the 

values determined from the shear modulus experiments . 23 ' 43 The magnitude of the relaxation time 

ri can be estimated crudely using the golden rule with r\ — J2e I (Ei\K\E 2 )\ 2 p(E)5(Ei — E 2 — 

E) where K is the lattice distortion and the density of states p(E) ~ E 2 / E 2 D . For E\ — E 2 = J, 

the exchange energy, and E D the Debye energy we find T\ju x ~ 4.10 -4 . In Fig. 10 we show the 

fit to the relaxation times T\, using T\obs = Tie + T x(T) for r = 8.3 x 1CT 9 s. and E = 1.8 

K. These values are to be compared with those inferred from the shear modulus experiments, 2123 

t = 2.3 x 1CT 9 s. and E = 1.7 K. 21 Although the values of r and E are a little different from 

those deduced by Beamish et al., 22 ' 23 the fit is remarkably good given the approximations that 

have been made. This interpretation of the results is consistent with the observations of Sasaki et 

al. 55 ' 56 who followed the temperature dependence of the NMR signals of 10 ppm 3 He in solid 4 He 

for densities such that the phase separation forms solid clusters of 3 He. While they were able to 

observe the solid clusters they were not able to detect the signal from isolated 3 He atoms above the 

phase transition at low temperatures, and they attributed this to a long spin-lattice relaxation time 

that lead to saturation and signal loss at the temperatures of interest for exploring the anomalies of 

the 4 He lattice. 



C. Effect of dislocations on relaxation times 

The most significant effect of the presence of dislocations on the NMR studies is that the 3 He 
atoms will pin the dislocation lines leaving some of the 3 He atoms located on the dislocation lines 
and the lines themselves essentially immobile at low temperatures. The number of 3 He atoms that 
are expected to be affected by this is very small. The localization of the 3 He atoms would lead 
to three distinct components of the NMR signals and their relaxation: (i) the contribution from 
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FIG. 10. (Color online) Comparison of the observed temperature dependence of the nuclear spin-lattice 
relaxation time T\ (for a sample with X3 = 16 ppm) with the dependence calculated for lattice relaxation 
due to the motion of the distortion of the lattice around impurities, using parameters inferred from the 
measurements of the shear modulus. 

atoms fixed on the dislocation lines, (ii) the interactions between fixed localized 3 He atoms with 
those tunneling freely in the solid, and (iii) the interactions between pairs where both are free to 
move in the lattice. The latter has been the focus of the discussions above. For those atoms that 
are essentially fixed the line shapes would be very broad leading to short values of T 2 and very 
long values for T\. The echoes would have more than one contribution and if the numbers of 
fixed atoms were significant the echo would have a central peak, corresponding to the spin-spin 
interactions between immobile atoms. We have not observed such a peak or a deviation from a 
single exponential decay for the dilute samples down to the lowest temperatures studied. 

This result is not surprising when one considers the number of dislocations present even in 
crystals of average quality. Paalanen et al. 51 and Iwasa et a/. 33 ' 58 show that the dislocation density 
A and the dislocation network length L N are related by AL 2 N = 0.2 and A = 10 6 cm~ 2 with 
L N ~ (3 — 6) x 10~ 4 cm. ~ 10 4 a . The concentration of 3 He atoms that would saturate the 
dislocation network is therefore estimated as x^ s ~ (j^) 2 ~ 10~ 8 . The dislocation density would 
need to be much larger than that expected for even average quality crystals 33 ' 58 if the pinned 3 He 
atoms were to be observable. 

In addition to dislocation lines, Balibar has pointed out 59 that grain boundaries are also likely 
to bind 3 He impurities - possibly affecting NMR relaxation times in a temperature-dependent 
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manner. A key parameter determining the importance of grain -boundary binding is the typical 
grain size Lq. We have no direct information on Lq for our blocked-capillary samples, although, 
for example, L G fa 0.5 mm was observed by Franck et al. 60 in annealed helium films. 

The grain size L G controls any possible effect of grain-boundary binding on the NMR response 
in two distinct ways: (a) the maximum concentration of 3 He that can be bound at grain boundaries 
is of order a /Lc, assuming full monolayer coverage with bound 3 He, and (b) the time scale 
for 3 He atoms to freely diffuse on and off grain boundaries is of order L G /D, where D is the 
diffusion coefficient. For Lq = 0.5 mm we calculate a /Lc = 0.7 ppm, much less than the 3 He 
concentrations used in the present experiments. For a diffusion constant D = 3.10 -6 cm 2 /s. 50 we 
find Lq/D falO 3 s., much shorter than the relaxation times T\ that we observe. The latter implies 
that the collisions between the 3 He impurities in the bulk sample dominates the relaxation process. 
Therefore, it appears that the grain size in our samples would need to be much smaller than was 
observed in films in Ref.[59] for the binding of 3 He atoms at grain boundaries to be the cause of 
the Ti anomalies we observe. 

D. NMR Relaxation for Droplets 

For all samples studied (except the 24 ppm sample) a sharp phase transition is observed at low 
temperatures below which the NMR amplitudes are independent of temperature (see Fig. 4) as ex- 
pected for the well-known phase separation into Fermi liquid droplets. 61 The observed phase sep- 
aration temperatures are in good agreement with the values calculated by Edwards and Balibar. 61 
The rate of formation of the droplets is very slow (typically 2-10 hours) 62 and great care must be 
taken to ensure that equilibrium is reached before measuring amplitudes and relaxation times near 
the phase separation. After phase separation, the relaxation times are observed to become temper- 
ature independent. 62 ' 63 The observed relaxation times are consistent with a relaxation that occurs 
at the wall of degenerate Fermi liquid droplets. Huan et al. 62 showed that the relaxation time was 
given by T\ = (N core /N wa u)Tx where N core and N wa u are the number of atoms in the core of the 
droplet and in the wall respectively. Tx is the intrinsic relaxation time at the wall due to the 3 He 
tunneling is given by tx = J3/M2 where J3 is the tunneling rate at the wall and M 2 is the NMR 
second moment. This estimate leads to values of T\ ~ 170 s. in good qualitative agreement with 
the observations (see Fig. 5). 
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V. CONCLUSION 



Measurements of the nuclear spin-lattice relaxation times for very dilute 3 He concentrations 
have shown the existence of pronounced peaks in the relaxation times at the same temperatures as 
those for which anomalies are observed in torsional oscillator and shear modulus measurements 
of solid 4 He. Less well-defined variations are observed for the nuclear spin-spin relaxation times. 
The detailed temperature dependences do not fit a model in which there is a well-defined phase 
transition to a supersolid or superfluid state where the critical fluctuations can induce dramatic 
changes in the relaxations times with critical exponents for the temperature dependence near the 
transition temperature. The observations are best described by the introduction of an additional re- 
laxation process in series with the usual tunneling-relaxation process and caused by the response 
of the lattice to the motion of the lattice distortions around the tunneling impurity atoms. The 
characteristic parameters for this model, the tunneling rate and the lattice excitation energy are re- 
markably close to those values deduced from measurements of the shear modulus. Further studies 
are needed at lower 3 He concentrations and lower magnetic fields to to create a better separation 
between the phase separation temperature and the temperatures for which the NMR and other 
anomalies are observed. 
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